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ABSTRACT 

During a galaxy merger, the supermassive black hole (SMBH) in each galaxy is thought to sink to 
the center of the potential and form a supermassive black hole binary; this binary can eject stars via 3- 
body scattering, bringing the SMBHs ever closer. In a static spherical galaxy model, the binary stalls 
at a separation of about a parsec after ejecting all the stars in its loss cone - this is the well-known 
final parsec problem. However it has been shown that SMBH binaries in non-spherical galactic nuclei 
harden at a nearly constant rate until reaching the gravitational wave regime. Here we use a suite of 
direct ./V-body simulations to follow SMBH binary evolution in both corotating and counterrotating 
flattened galaxy models. For N > 500A', we find that the evolution of the SMBH binary is convergent, 
and is independent of the particle number. Rotation in general increases the hardening rate of SMBH 
binaries even more effectively than galaxy geometry alone. SMBH binary hardening rates are similar 
for co- and counterrotating galaxies. In the corotating case, the center of mass of SMBH binary 
settles into an orbit that is in a corotation resonance with the background rotating model, and the 
coalescence time is roughly few hundred Myr faster than a non-rotating flattened model. We find 
that counterrotation drives SMBHs to coalesce on a nearly radial orbit promptly after forming a hard 
binary. We discuss the implications for gravitational wave astronomy, hypervelocity star production, 
and the effect on the structure of the host galaxy. 

Subject headings: Stellar dynamics - black hole physics - Galaxies: kinematics and dynamics - Galaxy: 
center. 


1. INTRODUCTION 

Supermassive black holes (SMBH)s, with masses be¬ 
tween 10 6 and IQ 1 0 solar masses, l i e at the h eart of nearly 
every galaxy (e.g.|Kormendy & Richstone|1995). With a 
few notable exceptions, these SMBHs dwell within stellar 
spheroids - spiral bulges, ellipticals, and SOs - and obser- 


erties of its host spherioid (e.g. Gebhardt et al. !2000| 

Fer- 

rarese & Merritt ;2000; McConnell & IVla 

2013 Gra 

iam 

& Scott|2015p. This seems to show that t 
the SMBH and its host are deeply tied, an 
observational and theoretical studies bear 

re evolution of 
d innumerable 
this out (ISher- 


man_eWiLJ[20M JVolonteri&^ellovary 2012; Hopkins & 

|Quataert|2010| |Micic et al. 120071 IDi Matteo et al.||2005|). 

'There is an emerging picture that SiVIBH hosts ail ro- 

tate to some degree. In terms of spiral bulges, rota¬ 
tional support is common. The early-type stars in our 
own gal actic center have a net counterrotat ion of 120 
km/sec (Genzel et al.|1996 Schodel et al. 2009). The pro¬ 
totypical pseudobulges are rapi dly-ro tating, aisky struc¬ 


tures (Kormendy & Kennic utt|[20 04) , an d most classical 


bulges show some rotation as well ( Gadotti 2012). 


disky and rotationally-supported (e.g. 

Faber et al. |1997 

Bender et al. 

1994 Cappellari et al. 

2007;). However, 


volume-limited census of the nearest 250 early-type 
galaxies ATLAS 30 fo und that an astonishin g 86 per¬ 
cent are fast rotators (Emsellem et al. 2011), with less 
than 3 percent of the sample being described as non¬ 
rotating. One of these non -rotators was thought to be 


the giant cD galaxy M87 (Emsellem et al. 2004 Kra- 


jnovic et al. 2011) 
this 


and yet new IFU data reveal that 
canonical non-rotator has an unmistakable 


kinematically-distinct core (Arnold et al. 2014 Emsellem 


even 



structure is quite complex. 

It is easy, theoretically, to expect that rotation is prac¬ 
tically ubiquitous in stellar spheroids. Gas-rich forrna- 


major mergers (Bender et al. 1992 

Khochfar & Silk 2009 

Bois et al. 

2011 

Tsatsi et al. 2015 

). indeed, a non- 


may be created excl usively 


minor mergers (Naab et al. 
clearly rotate, and since the 


throu gh a slew of gas-poor 


2014). Since SMBH hosts 
ink between the SMBH and 


its host is so well-established, any study of the dynamics 
with a SMBH host should consider rotation. This paper 
explores the evolution of binary SMBHs in a rotating, 
flattened stellar system, and its effect on the structure of 
the stellar host. 

Binary SMBHs are expected to form within the galaxy 
core after a merger, and if the two SMBHs coalesce, they 
are arguably the most powerful so urces of gravitational 
radiation in the Universe (Hughes 2003). One problem 
that plagues SMBHs is how they merge together. The 
binary can eject stars via 3-body scattering, bringing 
the SMBHs ever closer. Once the ejected stars extract 
enough energy from the binary orbit to shrink the sep¬ 
aration to roughly milli-parsec scales, gravitational ra- 
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diation dom inates, and the SMBHs coalesce (Begelman 


et al. 1980). However, the problem is that analytical 


calculations and simulations of static, spherical galaxies 
show that the binary’s orbital separation stalls bef o re the 


SMB Hs can plunge toward merger (see (Milosavljevic & 
Merri~tt]|2003) for a review). The root cause of this hang 


up is"simply a lack of low-angular momentum stars capa¬ 
ble of interacting with the binary via 3-body scattering. 
This theoretical bottleneck has become known as the in¬ 
famous “final parsec problem”. Recent work by several 


more realistically-shaped galaxy models (Berczik et al. 

2006: Khan et al. 

2011 

Preto et al.l 2011; Guaiandris 

& Merritt 2012 j K 

han et al. 2013p, because the stellar 


replenish t he SM BH b inary l oss cone (e.g. [Yu| 2002| (Mer¬ 


ritt & Poon 2004; Holley-Bockelmann & Sigurdsson 2006; 

Vasiiiev et al. ||2014 Li et al.||2014|). 

SMBH binary dynamics has been studied in spherical 
rotating systems, but the resolution of the simulations 
(fewer than 100000 particles) could not accurately track 
the evolution of the binary to the gravitational wave 
regime. Experiments suggest that we need of the or¬ 
der of 1 million particles to resolve the true evolution of 
the binary beyond the hard binary stage. These previous 
studies primarily looked at the evolution of binar y ecc en¬ 
tricity and inclination with respect to the host (Sesana 


et al.|20lT Amaro-Seoane et al.|2010 Wang et al.|2()l l ~ 
and discovered that the eccentricity of the binary in- 
creased dramatically in counterrotating systems, while 
in corotating systems, the binary tends to circularize. 

Here we study SMBH binary dynamics and evolution 
in a flat rotating galaxy model with b/a = 0 . 8 , which is 
quite a bit less flattened than a realistic galactic nucleus. 
For this model, we found an IV-independent SMBH bi¬ 
nary evolution in Khan et al. (2013). In this study we 
explore the dependence of SMBH binary evolution on 
particle number for various amounts of rotation of sur¬ 
rounding cusp. We also study the energy, angular mo¬ 
mentum and eccentricity evolution of the SMBH binary 
and estimate the coalescence times by scaling our model 
to various observed nearby galaxies. 

The paper is organized in the following way: In Sec¬ 
tion 2 we describe the initial conditions and numerical 
methods for our direct IV-body simulations. The results 
for the evolution of the SMBH binary in rotating axisym- 
rnetric galaxy models are explained in Section 3. Finally, 
Section 4 concludes, discussing caveats and future work. 


2. INITIAL CONDITIONS AND NUMERICAL METHODS 
2 . 1 . The host galaxies and their SMBHs 
To isolate the effect of rotation, we used the Flat 8 


model in Khan et al. (2013) (hereafter KH13) to generate 
our rotating models; this allows us to directly compare 
the SMBH binary orbital evolution results here to those 
in the identical non-rotating model. For this flattening 
ratio, we found that the SMBH binary evolved into the 
gravitational wave regime, and the evolution did not de¬ 
pend on particle number for N > 500 1\. 

We introduced rotation by flipping the z component of 
the angular momentum, L z , in the positive direction for a 
subset of particles which have negative L z . Our fiducial 
model, the A series, flips every particle with negative 



log r 




Fig. 1.— Evolution of our flat rotating model (A0) in isolation. 
Top Panel: Stellar density profile at various times. 7 = 1.0 is 
the reference theoretical profile, shown as a dashed line. Middle 
Panel: Evolution of intermediate to major (b/a) and minor to 
major (c/a) axis ratio as measured at the half mass radius. Bottom 
Panel: Ratio of rotational velocity to the 3-d velocity dispersion as 
a function of distance at various times. In isolation, the rotational 
support in the center decreases over time as the central SMBH 
increases the velocity dispersion. 

L-, while our B series only flips 50% of the negative 
L z particles, meaning that 75% of the particles have a 
positive L z . 

Our galaxy models have an inner density slope, 7 , of 
1.0 and minor to major axis ratio of 0 . 8 , measured at 
the half-mass radius. In model units, the total galaxy 
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TABLE 1 

Rotating Axisymmetric SMBH Binary Parameters 


Run 

N 

7 

c/a 

Rotation 

so 

1500 k 

1.0 

1.0 

none 

SI 

1000k 

1.0 

1.0 

none 

S 2 

500 k 

1.0 

1.0 

none 

A 0 

1500 k 

1.0 

0.8 

corotating 

Al 

1000k 

1.0 

0.8 

corotating 

A 2 

800 k 

1.0 

0.8 

corotating 

A 3 

500 k 

1.0 

0.8 

corotating 

A 4 

250 k 

1.0 

0.8 

corotating 

B 0 

1500 k 

1.0 

0.8 

0 . 75 % corotating 

B 1 

1000k 

1.0 

0.8 

0 . 75 % corotating 

B 2 

800 k 

1.0 

0.8 

0 . 75 % corotating 

B 3 

500 k 

1.0 

0.8 

0 . 75 % corotating 

B 4 

250 k 

1.0 

0.8 

0 . 75 % corotating 

CO 

1500 k 

1.0 

0.8 

counterrotating 

Cl 

1000k 

1.0 

0.8 

counterrotating 

C 2 

800 k 

1.0 

0.8 

counterrotating 

C 3 

500 k 

1.0 

0.8 

counterrotating 


Note. — Column 1 : Galaxy model. Column 2 : Number of 
particles. Column 3 : Central density slope 7. Column 4 : Axes 
ratio. Column 5 : Sense of galaxy rotation with respect to the 
SMBH binary initial orbit. 


mass is 1 and the SMBH at the center of each model 
has a mass of 0.005. Because w e adiabatically squeeze 


the mod el to generate its shape (Holley-Bockelmann et 
al.|20011 , the system can and does change its shape with 
radius as it adiabatically adjusts to a live and changing 
SMBH-embedded potential. In addition, the scale radius 
slowly shrinks from 1.0 to 0.5; we could resize our model 
such that the scale radius is 1.0 again, but we instead 
choose to adopt physical units to reflect this increased 
central density. Discussion on scaling of our model to 
observed galaxies is given in section [4] and Tabic [2] scales 
the model to physical units. 

To determine the stability of our rotating models, we 
ran them in isolation for 40 time units; this duration is 
about half the maximum evolution time of the SMBH 
binaries in this study. Figure [l] shows that both the den¬ 
sity profile and half-mass axes ratio remain very stable 
for whole duration of the run. The kinematics are also 
fairly stable; inside the radius of influence, however, the 
system does become hotter, with v/o decreasing from 
0.6 to 0.4 as the SMBH re-establishes its characteristic 
cusp in velocity dispersion. This decrease in rotational 
support affects only the innermost ~ 5000 particles, but 
since it is also in the region of the model that exhibits 
a less flattened, but triaxial shape, it may well be that 
this region contains or bits that are less stable to rota¬ 
tion (Deibel et al.|2011). The orbit content of this model 
will be a subject for future study. 

To explore the SMBH binary evolution in flat rotating 
galaxies, we introduce an equal mass secondary SMBH at 
a distance of 0.5 with 70 % of the galaxy’s circular veloc¬ 
ity at that initial separation. We investigated SMBH bi¬ 
nary orbits that are corotating and counterrotating with 
the sense of the galaxy rotation. Table [l] describes pa¬ 
rameters of our SMBH binary study. Note that we also 
evolved the SMBH binary in non-rotating spherical and 
flattened galaxy models with the same density profile to 
facilitate the comparison. 

The numerical methods and hardware used for this 
work is described in section 2.2 of KH13. 


3 . SMBH BINARY EVOLUTION IN FLAT ROTATING 
GALAXY MODELS 

Here we discuss the results of our numerical studies of 
SMBH binary evolution in flat rotating galaxy models. 
The top panel of Figure [2l shows the evolution of inverse 
semi-major axis for the models. We see that for N 
greater than 500K, the invers e semi major axi s evolution 
is independent of N, unlike in Vasiliev et al. (2014). To 
be conservative, we approach IV as high as 1.5 million, 
never used before in such a study. For reference, we also 
plot the 1 /a evolution of a flat non-rotating model with 1 
million from our previous study (KH13), as well as a 1.5 
million particle run in spherical galaxy (SO) model with 
the same density profile as our rotating galaxy model. 
We can see that in rotating flat models, the SMBH binary 
evolves at a rate considerably faster than in mere flat 
galaxy models. We see TV-independent evolution of the 
SMBH binary in flat rotating galaxy models for N as 
large as 1.5 million. This points to a potential stellar 
dynamical solution to the final parsec problem within 
flattened galaxy models. 

We also studied SMBH binary evolution in a galaxy 
with less dramatic bulk rotation; in this case, only 75 
% particles corotate with the massive binary (models B 
in table |T]) . The middle panel of figure [ 2 ] again shows 
that the model experiences ./V-independent evolution of 
1/a. for N > 8001G However the binary coalesces at 
a slower pace when compared to models A for the same 
particle number. For example, the convergent A models 
pass l/a= 1500 at 40 time units, while the convergent B 
models take until nearly 50 time units to pass this same 
point. 

The bottom panel of figure [2] shows the evolution in 
the semi-major axis for the SMBH binary in a counter¬ 
rotating orbital orientation (models C). Here, the binary 
orbit shrinks at slightly different rates than the previous 
models at various time intervals, though not with a clear 
trend. Af around 50 time units, the slope of 1/a line 
seems very similar for N > 500AT, and again we notice 
rapid evolution of 1/a when compared to flat and spheri¬ 
cal models. However we have doubts that we are captur¬ 
ing the SMBH binary evolution accurately, because there 
is a clear dependence on particle number at later stages. 
In section[4]we discuss why we believe that SMBH binary 
coalescence is achieved in our counterrotating models, 
despite the lack of convergence in the model suite. 

For consistency between runs, we calculate the hard¬ 
ening rates s for all our runs by fitting a straight line 
to the inverse semi major axis s = // (-) in the interval 
50-70 time units (fig [3]). Both corotating and counter¬ 
rotating binaries (runs A & C) have hardening rates of 
about s ~ 28 which is about 30% higher than mere flat 
models. The B runs have slightly lower values of s ~ 25 
when compared to (runs A & C). Overall, we find that s 
in both co- and counterrotating models is approximately 
4 times higher than in spherical models with exactly the 
same density profile, for our best resolved runs with 1.5 
million particles. 

For a spherical, homogeneous, isothermal background, 
the hardening parameter H is related to hardening rate 
s thr ough H = scr/Gp - for the full los s cone regime, H 
« 15 ( Quinlan|1996 Sesana et al.|2006 ). We calculate H 
by substituting values of p and a at the influence radius, 
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time 



time 



time 


Fig. 2. — Evolution of the semimajor axis of the SMBH binary 
for the A (top) B (middle) and C (bottom) model suites (see table 
fl}. For comparison, the grey lines represent non-rotating flattened 
(uppermost grey) and non-rotating spherical (bottommost grey) 
one million particle models. The opacity of the purple lines scale 
with the particle number - the most transparent line has the fewest 
particles in the suite. 


defined as sphere around the SMBH binary containing 
twice the mass of the binary in stars. Figure 0] shows 
that for runs A2, Al & AO with N > 800A;, where the 
evolution is independent of N, the value of H remains 
constant around 11. The value of H from these 7V-body 
simulations is within 70% of scattering experiments. For 
models B, again we see a constant value of H ~ 10 for 


30 


25 


20 


15 


10 


5 



Fig. 3 .— Hardening rates for all our numerical experiments. 
Grey points are for the spherical models (S); Purple represents 
models A; Green is for models B, and the red points are for models 
C. 
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Fig. 4. — Hardening parameter H for all our numerical experi¬ 
ments. Colors are as in figure [ 3 ] 


runs with N > 800/c. For spherical run SO with greatest 
particle number N, H ~ 1.9 almost 8 times smaller than 
what is predicted for a full loss cone in scattering exper¬ 
iments. We would like to point out that in our models, 
the background profile is not at all isothermal and there 
is also some ambiguity for where one should measure a 
and p to make a fair comparison with scattering experi¬ 
ments. With this in mind, it is very encouraging that H 
obtained from our study is well within a factor of two of 
idealized scattering experiments. 

Figure [5] shows the eccentricity evolution of SMBH bi¬ 
naries for all runs with N > 1 million. For corotat¬ 
ing SMBH binaries, the eccentricities are consistently 
small (e ~ 0.2). For counter-rotating SMBH binaries, 
the eccentricity approaches e ~ 1 as soon as the binary 
forms. This is consistent with the findings of|Sesana et 
|ah] (f20ll|). 

The dichotomy in eccentricity behavior is borne out in 
a difference in the evolution of the angular momentum 
loss. Figure [6] shows the angular momentum evolution of 
SMBH binaries for our best resolved co- and counterro¬ 
tating models. 

It is clear that in the counterrotating case, the angular 
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Fig. 5. — SMBH binaries eccentricity for models with N > 1 
million. SMBH binaries in co-rotating, flat and spherical models 
have very small values of eccentricity ~ 0.1 whereas eccentricity ap¬ 
proaches unity in counter-rotating models. Here, the transparency 
represents the particle number; opaque lines are the TV = 1.5 mil¬ 
lion runs, while the fainter lines show results for N = 1 million par¬ 
ticles. Orange lines show the counterrotating model; Green lines 
represent models B; Purple lines show models A, and the grey lines 
are for the spherical and non-rotating axisymmetric models. 



Fig. 7. — Motion of the center of mass of the SMBH binary in 
model AO (purple) and the non-rotating flattened model (grey). 
Fainter colors indicate earlier epochs. Note that for model AO, the 
binary center of mass settles into a roughly circular orbit about 
the galactic center with a radius of ~ 0.05 in model units, roughly 
the SMBH radius of influence. In contrast, the center of mass 
of the SMBH binary in the non-rotating model undergoes simple 
Brownian motion about the galactic center. 



time 



time 


Fig. 6. — Angular momentum evolution for corotating (top) and 
counterrotating (bottom) SMBH binaries. The green line is Lz, 
while the pink and brown lines are Lx and Ly, respectively. 



Fig. 8 .— Motion of the SMBH binary center of mass in model 
CO (purple) and the non-rotating flattened model (grey). The color 
transparency is as in figure [7] In our counterrotating model, the 
binary center of mass executes a random walk like the non-rotating 
model. 


momentum loss is much more rapid. As we see from fig¬ 
ure [2j the inverse semi-major axis evolution (and hence 
energy loss) is very similar for both co- and counterro¬ 
tating models. This faster loss of angular momentum 
translates into a rapid rise in eccentricity. 

We also investigated the center of mass motion of the 
binary in models A and C. Figure [7] shows the position 
of the center of mass throughout the run. 

The trajectory of the SMBH binary center of mass is 
strongly effected by rotation. For non-rotating flattened 
models, the binary center of mass exhibits a s mall ra n¬ 


dom w alk characteristic of Brownian motion (Chatterjee 
et al.|2002). On the other hand, the center of mass in the 


corotatmg system settles into a corotation resonance at 
the radius of influence, following a roughly circular orbit 
of radius nearly in the x-y plane. 

The counterrotating case, on the other hand, shows no 
binary orbital coupling, and the SMBH binary center of 
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TABLE 2 

Physical Scaling of our Models 


Galaxy 

M.(M 0 ) 

Th(pc) 

T(Myr) 

L(kpc) 

M(Mq) 

M87 

3.6 x 10 9 

460 

3.92 

3.68 

7.2 x 10 11 

NGC4472 

5.94 x 10 8 

130 

1.45 

1.04 

1.2 x 10 11 

N4486A 

1.3 x 10 7 

31 

1.14 

0.25 

2.6 x 10 9 


Note. — Columns from left to right; (1) Reference galaxy, (2) 
Observed SMBH mass, (3) Observed SMBH radius of influence, 
(4) time unit, (5) length unit, (6) mass unit 


mass undergoes a random walk very similar to the non¬ 
rotating case. 


4. SMBH BINARY COALESCENCE 

We also estimate the SMBH binary evolution for each 
case after the end of our simulation. We choose three 
Virgo cluster galaxies as reference to physically scale our 
models: M87, NGC4472, NGC4486A. In each case, the 
mass scale is set by the observed SMBH mass. For length 
scale, we set the influence radius of the SMBH binary- 
embedded galaxy model to the size of the observed in¬ 
fluence radius of reference galaxy model. NGC4472 may 
best represent our model density profile with its mod¬ 
erate central cusp, while the central core in M87 and 
the steep cusp of NGC4486A span the range of typical 
density cusps. Table [2] shows useful quantities for the 
physical scales in our models. 

Our technique for extrapolating the evolution of the 
SMBH binary beyond the endpoint of the simulatio n is 
explained in detail in section 4.3 of Khan et al. (2012b). 
We choose the runs with highest particle number in each 
model (AO, BO, CO) for this extrapolation technique, and 
the evolution is shown in figure [9] The top panel shows 
the SMBH binary evolution scaled to M87. 

In the non-rotating case, the SMBH binary coalesces 
in roughly 1.5 Gyr for this physical scaling. For corotat¬ 
ing models AO and BO, coalescence times are roughly 1.3 
and 1.1 Gyr. For the counterrotating run CO, we evolve 
the SMBH binary from 1/a = 1000 (see figure [2j) at a 
system time T = 27, and we only consider hardening by 
gravitational waves; this is because we are not certain 
that the scattering results converge when the binary or¬ 
bit is smaller than this. The SMBH binary coalesces in a 
mere 100 Myr - essentially immediately - due to its near 
radial eccentricity. The middle panel of figure [9] shows 
the SMBH binary evolution for NGC4472. Here, the 
SMBH binaries coalesce approximately two times faster 
than in M87; a case in point: the corotating model A0, 
the SMBH binary coalesces in roughly 500 Myr. Finally, 
the bottom panel of [9] scales to NGC4486A, and in this 
case the SMBHs merge in model A0 in about 1.5 Gyr 
while the SMBH binaries coalesce in almost 2 Gyr within 
the non-rotating model. Neglecting stellar hardening, 
the SMBH binary in the counterrotating case coalesces 
in roughly 2 Gyr. However, if we assume that SMBH 
binary reaches an asymptotic hardening rate in CO, coa¬ 
lescence happens immediately after a hard binary forms. 
Clearly, SMBH binaries coalesce faster in rotating flat¬ 
tened models, but the mechanism behind the coalescence 
is very different depending on the sense of rotation. In 
corotating models, the rapid coalescence is due to higher 
hardening rates, but for counterrotating models high ec¬ 
centricity drives the merger. Out of three representative 
galaxies, the SMBH binary coalescence time is shortest 





time (Myr) 


Fig. 9.— Complete orbital evolution of the SMBH binary from 
formation to coalescence for runs A0 (yellow), B0 (orange), CO 
(blue) and Flat (gray) when our galaxy model is scaled to M87 (top 
panel), NGC4472 (middle panel) and NGC4486A (bottom panel). 
In the bottom panel, the very different counterrotating timescales 
either include (short timescale with fainter hue) or exclude (longer 
timescale with stronger hue) stellar hardening. 

for NGC4472. 

We show the coalescence time for different runs and 
physical scaling in Table |3j Technically, the clock starts 
here when two SMBHs form a pair with a separation 
~ 10 influence radii (presumably after a major galaxy 
merger), and ends with the coalescence of the binary from 
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TABLE 3 

Physical Scaling of our Models 


Run 

S 

H 

e 

T7.M87 

77,4472 

77,4486 A 

A0 

27.7 

10.87 

0.17 

1.17 

0.57 

1.52 

B0 

23.90 

9.20 

0.06 

1.27 

0.62 

1.67 

CO 

23.6 

8.66 

0.99 

0.11 

0.04 

2.04(0.14) 


Note. — Columns from left to right; (1) SMBH binary evolution 
run, (2) SMBH binary hardening rate, (3) Hardening parameter H, 
(4) SMBH binary eccentricity. (5) Coalescence time (in Gyr) when 
our model is scaled to M87, (6) NGC4472 (7) and NGC4486A. 



The 3-body scattering phase removes roughly 1.3 times 
the binary SMBH mass, scouring the density cusp out to 
about 1.1 kiloparsec if scaled to the M87 core. 

We found that when the SMBH binary and the galaxy 
are corotating, the eccentricity remains low at approx¬ 
imately 0.1, while counter-rotation acts to pump the 
SMBH binary eccentricity ud to nearly one during the 



Fig. 11.— This is the initial (orange) and final (red) 3- 
dimensional velocity dispersion profile of model A, showing a sharp 
spike in the velocity dispersion at the end of the run. 


Fig. 10.— This is the initial (green) and final (orange) density 
profile of model A, showing the clear mass deficit out to r ~ 0.3 of 
1.3 times the binary SMBH mass. 

gravitational radiation. 


5. STRUCTURE OF THE MERGER REMNANT 

In this section, we go over the imprint of the binary 
black hole merger on the galaxy remnant. As expected, 
figure [lO] shows that the density cusp is scoured out 
by 3-bod y scattering as the black holes coalesce (Gra¬ 
ham 2004). One puzzling result can be seen in figure 
ECO where the final velocity dispersion spikes; this is in 
contradiction to the dip in veloci ty dispersion expected 


from the stellar hardening phase (Meiron & Laor 


2013), 


and fully consistent with the velocity dispersion cusp ex¬ 
pected in an equilibrium SMBH-embedded nucleus. Fur¬ 
ther study is needed, using simulations with shorter snap¬ 
shot output cadence, to help pinpoint the occurrence and 
longevity of this potential kinematic signature of 3-body 
scattering. 


6 . CONCLUSION 

We investigated the effect of bulk rotation on SMBH 
binary coalescence in N-body-generated flattened galaxy 
models. Overall, we found that rotation drives the 
SMBH binary more efficiently through the 3-body scat¬ 
tering phase, resulting in coalescence timescales that are 
between 3 to 30 times faster than the same non-rotating 
model for co- and counterrotating models, respectively. 


inspiral phase. Such a high eccentricity enhances the coa¬ 
lescence of t he SMBH binary, as is seen in man y previous 
studies (e.g. Sesana 2010 Khan et al. 2012a). Though 


we caution that the eccentricity behavior is not conver¬ 
gent even for 1.5 million particles, we suspect that the 
eccentricity will remain high in the convergent regime; 
when the SMBH is counterrotating, the abundance of 
retrograde orbits can extract angular momentum from 
the bina ry very efficiently, and se cular dynamical anti¬ 
friction (Madigan & Levin 2012) torques the orbit so 
that it bleeds angular momentum! A systematic study is 
needed to gauge the degree of counterrotation and binary 
orbital plane alignment needed to pump the eccentricity 
into the nearly radial regime; if the binary eccentricity is 
very sensitive to minor degrees of counterrotation, then 
few mergers will linger in the 3-body scattering stage. 
For such high eccentricities, we should expect residual ec¬ 
centricity to persist into the last few orbits in the gravita¬ 
tional wave regime; this will have profound implications 
for gravitational wave detection using waveform template 
matching. 
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